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ABSTRACT 

We build, as far as theory will permit, self consistent model H II regions around central clusters 
of aging stars. These produce strong emission line diagnostics applicable to either individual H II 
regions in galaxies, or to the integrated emission line spectra of disk or starburst galaxies. The models 
assume that the expansion and internal pressure of individual H II regions is driven by the net input of 
mechanical energy from the central cluster, be it through winds or supernova events. This eliminates 
the ionization parameter as a free variable, replacing it with a parameter which depends on the ratio 
of the cluster mass to the pressure in the surrounding interstellar medium. These models explain why 
H II regions with low abundances have high excitation, and demonstrate that at least part of the 
warm ionized medium is the result of overlapping faint, old, large, and low pressure H II regions. We 
present line ratios (at both optical and IR wavelengths) which provide reliable abundance diagnostics 
for both single H II regions or for integrated galaxy spectra, and we find a number that can be used 
to estimate the mean age of the cluster stars exciting individual H II regions. 

Subject headings: galaxies: general - galaxies: star formation rates - galaxies: abundances - galaxies: 
starburst -ISM: H II -ISM: abundances- 



1. INTRODUCTION 

Much of what we have learnt about the chemical evolu- 
tion of the Universe, or of individual galaxies throughout 
cosmic time has been gleaned from the study of the spa- 
tially unresolved emission line specta of distant galaxies 
(Steidel et al. 1996; Kobulnicky, Kennicutt & Pizagno 
1999; Kobulnicky & Zaritsky 1999; Kewley & Kobul- 
nicky 2005). This emission line spectrum arises from 
the ensemble average of H II regions within the galaxy, 
characterized by a wide range of physical parameters. 

For example, in disk galaxies, the chemical abundance 
of heavy elements falls continuously from centre to edge. 
In like manner, both the pressure and density of the in- 
terstellar medium in the disk drops exponentially with 
radius. For a given input of mechanical energy by the 
central star cluster, this means that the outermost H II 
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regions expand more rapidly, and are both larger and 
have lower internal pressure at a given age. However, the 
properties of the central cluster of a given mass change 
systematically with chemical abundance. Massive stars 
with higher chemical abundances lose more mass prior to 
supernova explosion, have stronger stellar winds, evolve 
more rapidly, and have lower effective temperatures while 
they are on the Main Sequence. 

Finally, the temporal evolution of individual H II re- 
gions is extremely important in determining the inte- 
grated emission line spectrum. Leitherer, Gruenwald & 
Schmutz (1992) demonstrated the importance of aging of 
the OB stars in lowering the number of ionizing photons 
as the most massive stars evolve away from the Zero Age 
Main Sequence (ZAMS) to become supergiants. Later, 
these stars enjoy a brief resurgence in their ionizing pho- 
ton and mechanical energy production as they become 
Wolf-Rayet stars. 
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With so many temporal and physical variables, it 
is perhaps not surprising that the nebular modelling 
community (ourselves included) have in the past been 
reduced to making severe, even gross, oversimplifi- 
cations in attempts to model the integrated strong 
emission line spectra of galaxies. Indeed, frequently 
only a single line ratio has been used to determine 
the metallicity of extragalactic H II regions or even 
of whole galaxies. This is the famous R 2 3 ratio; 
([O II] A3727,29+[0 III] AA4959,5007)/H/3. 

This ratio was first proposed by Pagel et al. (1979). 
The logic for its use is impeccable, since it uses the two 
strongest lines of the strongest coolant of H II regions, 
and it should therefore be sensitive to the total oxygen 
abundance. By using two stages of ionization, it ac- 
counts for the emission in the bulk of the ionized neb- 
ular volume. However, the calibration of this ratio in 
terms of the abundance has proved to be very difficult 
and many different and often contradictory calibrations 
of R23 are available, including Pagel et al. (1979); Pagel, 
Edmunds & Smith (1980); Edmunds & Pagel (1984); Mc- 
Call, Rybski & Shields (85); Dopita & Evans (1986); 
Torres- Peimbert, Peimbert & Fierro (1989); Skillman, 
Kennicutt & Hodge (1989); McGaugh (1991); Zaritsky, 
Kennicutt & Huchra (1994); Pilyugin (2000); Kewley et 
al. (2002b). 

The basic cause of this difficulty in calibrating the 
R23 ratio is common to many ratios formed from the 
flux of optical forbidden line divided by the flux of a 
recombination line of hydrogen. Initially, as Z is in- 
creased, the forbidden line increases in flux, and the line 
ratio increases. However, as the abundance increases 
further, the cooling by forbidden lines lowers the elec- 
tron temperature making it more difficult to collision- 
ally excite the optical forbidden lines. Eventually there 
comes a point at which an increase in abundance is 
matched by a decrease in collisional excitation and the 
line ratio reaches its maximum. Any further increase 
in Z then leads to a decrease in the relative strength 
of the forbidden line. R23, like other optical diagnos- 
tics is therefore a two- valued function of Z. Because 
infrared lines have lower thresholds for collisional excita- 
tion, they remain monatonic functions of Z up to much 
higher mctallicities. Recognizing this physics, the S23 ra- 
tio; ([S II] AA6717,31+[S III] AA9069,9532)/H/3, or even 
more complex ratios, have been used either in the place 
of or to supplement the R23 ratio (Dennefeld & Stasinska 
1983; Kennicutt & Garnett 1996; Oey et al. 2000). 

The calibration problems of the R23 ratio are made 
worse by the scatter, at fixed R23 and within each given 
branch, both between models and observational data, 
and amongst different models. The observational ma- 
terial is affected by uncertainty in the reddening correc- 
tion applying to the [O II] AA3727,29A , since classical 
extinction laws are usually used, rather than the attenua- 
tion laws which more correctly apply to extended objects 
(Fischera, Dopita & Sutherland 2003; Fischera & Dopita 
2005). This point will be discussed more fully below. 
As far as the theoretical models arc concerned, scatter 
in the relationship is easily generated, since the stellar 
atmospheres used for single stars can show very wide 
differences in the number of ionizing photons produced 
per unit frequency (Morisset et al. 2004). 

Moving beyond the use of single line ratios, a more so- 



phisticated analysis is made possible by the understand- 
ing that the strong line emission spectrum of an individ- 
ual H II region is controlled by three physical parameters 
(Dopita & Evans 1986). These parameters are, the ef- 
fective temperature of the exciting stars, the chemical 
abundance set (or "metallicity", Z) and the ionization 
parameter, defined either as the ratio of the mean ion- 
izing photon flux to the mean atom density, q or in its 
dimensionless form as the ratio of mean photon density 
to mean atom density U — q/c. Thus, for a given stellar 
input spectrum a U : Z grid of models should provide 
the ideal means of interpreting the observational data 
(Kewley et al. 2002b). 

This analysis leaves open the issue of the age distribu- 
tion of the exciting stars. The metallicity and the age of 
the exciting stars determines the energy distribution of 
the EUV photons. In the modelling hitherto, this prob- 
lem is usually treated in two limits; the instantaneous 
burst approximation (for which zero age is often, but 
not universally assumed) or the continuous star forma- 
tion approximation in which stars are assumed to be born 
continually within the ionized region (see Kewley et al. 
(2001), and references therein). All this modeling effec- 
tively assumes a single H II region of a given metallicity 
and ionization parameter. Whilst such an approach is 
acceptable in attempts to understand the emission line 
spectrum of individual H II regions, it is certainly inad- 
equate in attempts to interpret the integrated spectrum 
of a whole galaxy, or even of a complex of many H II 
regions. 

The advent of the Starburst 99 v. 5 (Leitherer et al. 
1999) stellar spectral synthesis code finally allowed the 
construction of fully self-consistent models for H II re- 
gions. For a given stellar initial mass function (IMF), 
this code delivers not only the pan-spectral distribution 
of photons, including ionizing photons as a function of 
age, but also provides tables of the time-dependent me- 
chanical energy input of the stars including OB-star stel- 
lar winds, red giant winds, Wolf-Rayet winds and the 
energy input from supernova explosions. 

In the first paper of this series (Dopita et al. (2005), 
hereinafter referred to as Paper I), we computed the (one 
dimensional) evolution of H II regions around clusters of 
a single mass. We demonstrated how the size : age rela- 
tionship is essentially determined by the density (or pres- 
sure) in the surrounding interstellar medium (ISM) . This 
in turn determines the dust 'temperature' described by 
the wavelength of the peak of the far-IR dust re-emission 
feature. 

In the second paper of this series, Dopita et al. (2006), 
in press, hereinafter referred to as Paper II), we inves- 
tigated the role of metallicity and of the cluster mass 
function in controlling the excitation of the populations 
of H II regions in galaxies, be they disk or starburst. 

We demonstrated that the ionization parameter of H II 
regions, previously treated as a free variable, is deter- 
mined at any time by the instantaneous ratio of the ion- 
izing photon flux to the mechanical energy flux of the 
central stars. This is because the HII region is evolving 
as a mass-loss bubble pressurized by the combined ram 
pressure of the stellar wind and supernova explosions, 
so HII region density and HII region radius are closely 
coupled. Although this conclusion was obtained using a 
simple 1-D spherical evolution model for the H II region, 
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it should remain valid in more complex geometries, since 
wherever un- ionized or molecular inclusions exist, these 
will be pressure confined by the hot shocked stellar wind 
gas, and this gas pressure is coupled with the overall size 
of the H II region, which itself determines the dilution of 
the radiation field of the central stars. 

Because the effective temperature of the exciting stars 
increases as metallicity is lowered, but the mechanical en- 
ergy flux in the stellar winds decreases towards to lower 
metallicity, low abundance HII regions are characterized 
by higher ionization parameters, a fact that had been fre- 
quently noted by observers, but which had not previously 
been explained in a satisfactory manner. 

Using these insights, in this third paper of the series, 
we proceed to the next stage of sophistication in model- 
ing - the construction of time-dependent photoionization 
models of individual H II regions, and the construction 
of ensemble averages of aging H II regions to provide 
galaxy-wide averaged spectra at a given metallicity. Such 
models applied to individual H II regions make it possi- 
ble to estimate the ages of the exciting stars from the 
positions of the observations on theoretical H II region 
isochrones. Previously, this could not be done, since the 
ionization parameter and the effective temperature of the 
cluster (determined by the ageing of its stars) can both 
produce similar effects on the emitted spectrum. This 
degeneracy is raised thanks to the self-consistent geome- 
try that we have developed, which includes a dynamical 
evolution of the H II region consistent with the properties 
of the stellar wind generated by the central cluster. 

In this paper, we apply these new techniques to pro- 
duce line ratio diagnostics which enable the independent 
determination of both the stellar ages and metallicitics 
of the exciting clusters for individual H II regions. These 
should prove useful in the analysis of abundance gra- 
dients in resolved disk galaxies. More importantly, we 
develop line ratio diagnostics for the determination of 
abundances in ensembles of H II regions in starburst and 
normal galaxies. This should greatly assist in analyses 
of the strong emission lines of the more distant and un- 
resolved galaxies, observed either at optical wavelengths 
or in the IR. 

2. MODELS 
2.1. Codes 

We have used the Starburst 99 v. 5 code in its latest 
(2005) version to compute the pan-spectral distribution 
of clusters of stars with a piece-wise fit to a Miller-Scalo 
initial mass function (IMF) Miller & Scalo (1979) be- 
tween 0.1 and 120 M Q . The parameters of this fit are 
the same as given in Miller & Scalo (1979) and were also 
tabulated by us in Paper II. Either the use of a Miller- 
Scalo IMF or of a Salpeter IMF would not noticeably 
affect the line ratio diagnostics presented in this paper, 
since the slopes of the IMF are almost identical above 
10 M Q . However, the choice of the IMF has a much 
larger effect on the total mass of the cluster, since this is 
determined by the choice of the power law below 10 M Q . 

Stellar atmosphere models for stars with plane paral- 
lel atmospheres are based on the Kurucz (1992) models 
as compiled by Lejeune (1997). The fully line-blanketed 
Wolf-Rayet atmosphere models of Hillier & Miller (1998) 
and the non-LTE O-star atmospheres of Pauldrach, Hoff- 
man, & Lennon (2001) have been incorporated into Star- 



burst 99 v. 5 as described in Smith, Norris & Crowthcr 
(2002). The flux distributions following from these at- 
mospheres are coupled to the stellar evolution models by 
the Geneva group. We used the so-called "high-mass- 
loss" tracks, which provide the best match to the ob- 
served stellar inventory in the Hcrtzsprung-Russcll dia- 
gram (Leitherer et al. 1999). 

For the Initial Mass Function (IMF) we used piece- wise 
power-law fits to the Miller & Scalo (1979) mass func- 
tion, although we note that their paper, Miller & Scalo 
(1979) represented their IMF as a truncated lognormal 
distribution (with its maximum at zero). However, they 
also represented it as a three-segment broken power law, 
which is the form of the IMF which Starburst 99 cur- 
rently accepts. The logarithmic slopes of the piecewise 
power-law fits, T, and the mass range applicable to each 
segment were given in Table 1 of Paper II. 

For a given cluster mass, we used the stellar wind and 
supernova power output, as tabulated by the Starburst 
99 v. 5 code, to solve for the radius of the H II regions as 
a function of time given the pressure or, equivalently, the 
mean density in the ISM. This requires a Runge-Kutta 
integration of the equation of motion of the swept-up 
shell, as was done and described in both Paper I and 
Paper II. 

The pressure in the H II region is the same as in the 
shocked stellar wind gas. This is given by the classical 
Castor, McCray (1975) theory; 

( 7 /J250_\ 4/15 / L mcch (t) \ 2/3 M ro g n 

U (38507r) 2 /5 ^3087ry \ fxm H n J r(t) 4 /3 

(1) 

where £ mC ch(i) is the instantaneous production of me- 
chanical energy by the cluster stars, r(t) is the radius, 
no is the mean atom density in the surrounding ISM 
and n is the mean atom density in the ionized gas and 
S(t) = n/riQ is the instantaneous ratio of the density in- 
terior to the bubble to the density in the surrounding 
ISM. 

In Paper II, we showed that the ionization parameter 
at the contact discontinuity in the swept-up shell of ISM, 
which can be treated as the inner boundary of the H II 
region depends largely on the instantaneous properties 
of the exciting cluster stars; 

q(t) ^5(t) 3 / 2 S*(t)/L mcch (t). (2) 
where S 1 * (t) is the instantaneous flux of ionizing photons 
from the central cluster. 

The presence of the S(t) factor in the above equation 
provides a weak coupling between the ionization param- 
eter and both the pressure in the ISM, Pq, and the mass 
of the central cluster, M c i. Together, these determine 
the strength of the outer shock of the mass-loss bubble 
and therefore the compression factor through it. The ap- 
propriate scaling factors are q cx P 1/5 and q oc M 1 / 5 . 

1 /5 

Thus, once Po and q oc MJ are fixed, the ionization 
parameter is determined. 

Using this dynamical model to derive the instanta- 
neous size and ionization parameter of the H II region, 
the latest version of our code Mappings Illr (Sutherland 
& Dopita 1993; Dopita et al. 2002; 'Groves, Dopita & 
Sutherland 2004) with the physical parameters defined 
in Paper I was used to compute photoionization mod- 
els. The resulting model isobaric dusty H II regions were 
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used to produce a table of emission lines, described in 

the Appendix, and given fully the Electronic version of TABLE 1 - Adopted solar chemical abundance ratios, log[n x /n H ], and 
the paper. These enable us to compute the line ratio DEPLETION FACT0RS ( >■ 
diagnostics described in this paper. 

2.2. Metallicities & Depletion Factors 

The solar abundance set adopted here is taken from 
Asplund, Grevesse & Sauval (2005). This incorporates 
the results of many papers which together have provided 
a self-consistent re-calibration of the solar abundance. 
In many cases, these abundances are nearly a factor of 
two lower than those used previously (Anders & Grevesse 
1989). The corresponding logarithmic abundance ratios 
of the elements with respect to Hydrogen are given in ta- 
ble 1. For the purposes of the modeling, when we refer to 
solar metallicity, Zq, we are using this set of abundances. 

We note here that as a consequence of the re-scaling 
of the solar abundances (which is not yet reflected in the 
stellar models), the Starburst 99 "solar" abundance set 
is inconsistent with the new Asplund, Grevesse & Sauval 
(2005) abundance set. The Starburst 99 code uses ab- 
solute metallicities of Z = 0.001,0.004,0.008,0.02 and 
0.04, which makes the Starburst 99 "solar" metallicity, 
Z = 0.02, higher than the Asplund, Grevesse & Sauval 
(2005) "solar metallicity" , Z = 0.016. Whilst this mat- 
ter very little in terms of the number of ionizing photons 
emitted by the stars, it will mean that the computed stel- 
lar UV photon field will slightly softer than it should be, 
given the nebular abundance set. 

The depletion factors from the gaseous phase deter- 
mine the composition of the dust in the models. For these 
we have took the measured gas-phase abundances in the 
local interstellar cloud by Kimura, Mann & Jessberger 
(2003), and used these to infer a set of local depletion 
factors consistent with the Asplund, Grevesse & Sauval 
(2005) abundance set. The adopted depletion set is also 
given in table 1. For models with non-solar abundances, 
the depletion factors were held constant at solar values. 
This is equivalent to taking the dust : gas ratio to be 
proportional to metallicity. 

Three elements do not scale simply with the metallic- 
ity. First, because of its high initial abundance as a result 
of nuclear burning in the Big Bang, Helium scales only 
weakly with metallicity. For this element, we assume a 
primary nucleosynthesis component in addition to its pri- 
mordial value. From Russell & Dopita (1992) and Pagel 
et al. (1992) we infer the empirical relationship: 
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Fig. 1. — The fit given by equation 5 is shown along with the HST 
observations of H II regions made by Garnett et al. (1999). The 
solar abundance from Asplund, Grevesse & Sauval (2005) is also 
marked with a symbol. The new solar abundance set eliminates 
the systematic difference between the solar and the H II region 
abundances which was noted by Garnett et al. 



form as for N: 

C/H = 6.0 x 10- 5 Z/Z Q + 2.0 x 10" 4 Z/Z Q i 
The fit to the observations is shown in figure 1. 
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2.3. 



He/H = 0.0737 + 0.024 x Z/Z . 



(3) 



For both Nitrogen and Carbon there is clear evidence 
that these are both a primary nucleosynthctic element, 
dominant at low metallicity, and a secondary nucleosyn- 
thetic element once higher abundances are reached. Al- 
ternatively, both are produced in part by dredge-up in 
intermediate mass stars, causing a delay in the onset of 
enrichment of these elements from this source. For Nitro- 
gen we have used a modified version of the form proposed 
by Groves, Dopita & Sutherland (2004); 



N/H = 1.1 x IQ- 5 Z/Z Q + 4.9 x IQ~ 5 Z/Z Q [ 



(4) 



For Carbon, we have used the HST data on the Carbon 
abundances of H II regions as compiled by Garnett et al. 
(1999). This provides a fit of the same mathematical 



Parameters of the Photoionization Models 

The Starburst 99 v. 5 code allows us to study line 
ratios at five (fixed) absolute metallicities of Z = 
0.001, 0.004, 0.008, 0.02 and 0.04. For each of these abun- 
dance sets we have computed several families of H II 
region models with ages 0.2, 0.5, and 1.0 Myr and in 
further 0.5 Myr steps up to a maximum age of 6.5 Myr. 
By this time, more than 96% of ionizing photons have 
been emitted, the H II regions have faded by factors of 
between 5 and 12 from their peak luminosity, and their 
surface brightnesses will have faded by much greater fac- 
tors than these. 

The ionization parameter is not a free variable in these 
models. The geometry of the gas with respect to the 
ionizing stars is instead determined by the equation of 
motion of the expanding H II region under the driving 
pressure of the stellar mass-loss and supernova explo- 
sions. This pressure also determines the density of the 



ionized plasma. As was shown in Paper II, the equation 
of motion and the pressure together lead us to infer that 
there is a weak coupling between the ionization param- 
eter and both the pressure in the ISM, P , and he mass 
of the central cluster, M c /. Together, these determine 
the strength of the outer shock of the mass-loss bubble 
and therefore the compression factor through it. The 

appropriate scaling factors are q oc P Q and q oc M*/ 5 . 

The ratio, TZ = (M cl /M )/(P o /fc), with P /k, mea- 
sured in cgs units (cm _3 K), uniquely determines the run 
of ionization parameter with time. It is this variable that 
replaces the ionization parameter and other geometrical 
considerations in our models. In figure 2 we show the run 
of the computed ionization parameter q as a function of 
time and ISM pressure for a solar abundance cluster with 
M cl = 3 x 1O 3 M and having a Miller-Scalo IMF. The 
(dimensional) ionization parameter q is related to the 
more frequently used (dimensionless) ionization parame- 
ter U by U — q/c. 

The run of q with time is computed using the Starburst 
99 output, and solving for the time dependent radius and 
internal pressure of the mass-loss bubble as described 
above. In the early phases (t < 2 Myr) the ratio of 
ionizing photon flux to mechanical energy flux is almost 
constant, and the changes in q mostly reflect the size evo- 
lution of the mass-loss bubble. After 2 Myr, the ionizing 
photon flux rapidly decreases as massive stars evolve to 
red supergiants, and then the mechanical energy flux in- 
creases as a result of Wolf-Rayet stellar winds. Both of 
these lead to a strong decline in q. At around 3.5 Myr, 
supernova explosions add to the internal pressure of the 
bubble, sharply decreasing the q in the ionized gas and 
re-accelerating the expansion of the mass-loss bubble. At 
later times, the progressive death of the high-mass stars 
leads to further declines in the ionizing photon flux, but 
the decline in q is slower because the internal pressure 
of the mass-loss bubble also declines, partially offsetting 
the decline in ionizing photon flux. 

The theoretical ionization parameter rather strongly 
dependent upon the chemical abundance. It is driven 
by two factors. First, at higher abundance, the stel- 
lar wind has a higher opacity and therefore absorbs a 
greater fraction of the ionizing photons, reducing the q 
in the surrounding H II region. Second, the atmosphere 
scatters the photons emitted from the photosphere more 
efficiently when the atmospheric abundances are higher, 
leading to a greater conversion efficiency from luminous 
energy flux to mechanical energy flux in the stellar wind 
base region. This also leads to a diminution of q in the 
surrounding H II region. These factors acting together 
provide a sensitivity to mctallicity of q oc Z~ 0,8 , approx- 
imately. 

As described above at any age t, the instantaneous 
ionization parameter - which determines the excitation 
of the H II region - is governed by the ratio of the mass 
of the cluster and the pressure in the surrounding ISM. 
Since cluster masses may vary between ~ IOOMq and 
10 6 Mq, and the likely range of ISM pressures are be- 
tween 10 4 < Po/k < 10 7 cm~ 3 K it follows that reason- 
able values for TZ fall within the range —6 < log TZ < +2. 
We have therefore computed, for each mctallicity, an 
age sequence of model H II regions having log TZ = 
—6, —4, —2, and +2. For a given metallicity, the in- 
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Fig. 2. — The computed run of ionization parameter at the inner 
edge of the H II region as a function of time and ISM pressure 
logPo/k (cm - 3 K) for a cluster having a mass of 3 X 10 3 Mq and 
assuming a Miller-Scalo IMF. A given choice of the variable 1Z = 
(M c i/ Mq) / (Po / k) defines a unique run of the ionization parameter 
with time. 



stantaneous value of the ionization parameter depends 
on TZ 1 / 5 , as we have shown in Paper II. Note that a par- 
ticular choice of 1Z is not tied to any particular cluster 
mass. Thus, a cluster of 10 5 M would be characterized 
by log TZ ~ 0.6 in the solar vicinity, and log 1Z ~ —1.0 in 
a high-pressure starburst environment. 

It is worth noting at this point that this parameteriza- 
tion of the problem of evolving H II regions only remains 
valid for as long as stochastic variation in the cluster 
membership may be ignored. This means that the ex- 
citing cluster has to be massive enough that the upper 
IMF is adequately populated with stars so as to avoid 
stochastic variations in either the mass to light ratio or 
the effective temperature of the cluster stars. This effec- 
tively limits the cluster mass for which these models are 
valid, M c i > 1O 3 M . 

Together, the two parameters Z/Z Q and log TZ define 
a unique H II region spectrum, with one caveat, namely, 
that if the density of the ionized gas becomes too high, 
then collisional de-excitations of forbidden lines will al- 
ter the emergent spectrum. This only becomes a con- 
cern at the very highest values of ISM pressure cither 
within ultra-compact H II regions (as computed by Do- 
pita et al. (2006a) or else in starburst environments with 
Po/k > 10' cm~ 3 K. For the current purposes, we have 
computed the particular case of Po/k = 10 6 cm~ 3 K, 
which constrains the electron density in the H II region 
below ~ 100cm -3 except at very earliest times and in 
the models of highest mctallicities, where it may reach 
~ 10 3 cm -3 . 

For each model, we have prepared tabular data for all 
the strong lines between the Lyman Limit and 88/im, 
from which we form H/3 flux weighted time averaged spec- 
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tra to represent the integrated spectrum for a given Z/Zq 
and log TZ. The set of strong spectral line intensities with 
respect to H/3 are available on-line. The form of the ta- 
bles is given in the Appendix, Tables 2-7. The code 
also returns the intensities of very many more fainter 
lines, which are not listed here. 

3. H II REGION EMISSION LINE DIAGNOSTICS 

3.1. Optical Lines 

3.1.1. Veilleux & Osterbrock Diagnostics 

The diagnostics that have been most frequently used 
in the classification of the nature of nebular excitation 
are those proposed by Veilleux & Osterbrock (1987) 
(hereinafter V&O). In these, the [N II] A6584A/Ha the 
[S II] AA6717,3lA/ Ha or the [O I]A6300/Ha arc plotted 
against the [O III] A5007A/H/? ratio. 

These V&O ratios have the virtue of only using lines 
which are close together in wavelength in forming the 
ratio, so errors due to uncertain reddening corrections 
are avoided. These plots also nicely separate H II regions 
from regions excited by active galactic nuclei (Seyfert or 
LINER nuclei) and from shock-excited objects. The H II 
regions form an extremely tight sequence, particularly in 
the diagram using the [N II] A6584A/Ha ratio. 

The results of our modelling for the [N II] A6584A/Ha 
vs. the [O III]/H/3 ratio is shown in figure 3. Be- 
cause of the difficulty of representing so many theoreti- 
cal models on the same figure, this is presented here in 
color, and in four panels. The first panel shows mod- 
els of all ages for each metallicity, color coded for clar- 
ity. The model series corresponding to log TZ = — 2 is 
drawn bold, since this is our best guess of the most likely 
value of log U in disk galaxies (10 2 < M cl /M Q < 10 4 and 
10 4 < P/k < 10 6 cm~ 3 K). The remaining three panels 
show the complete grid of models for log TZ = 0, —2 and 
-4, respectively. Isochrones are plotted for 0.2, 1.0, 2.0, 
3.0 and 4.0 Myr. 

In this figure, the observational datasets for individual 
H II regions are presented as crosses. The data are drawn 
from van Zee et al. (1998); Kennicutt & Garnett (1996); 
Dennefeld & Stasinska (1983); Walsh & Roy (1997); Roy 
& Walsh (1997) and the extension of the data set to low 
abundances is from Pagel et al. (1992). 

The reason why the observed H II regions on a 
[N II] /Ha vs. the [O III]/H/3 form such a narrow se- 
quence is quite evident. First, to be bright, and there- 
fore selected for observation, observed H II regions should 
have an age less than ~ 3 Myr. In this age range, the 
isochrones of all models of different abundances fall into 
the very narrow strip defined by the observations. Al- 
though it depends somewhat on the assumed value of 
log TZ, it is clear that no observed H II region has an age 
greater than 3 — 4 Myr, consistent with the presumption 
of strong observational selection against old H II regions. 

This figure should be compared with Figure 2 of Dopita 
et al. (2000), in which a grid of models in Z/Zq and q, 
the ionization parameter were presented. It is clear that, 
in order to have a low metallicity H II region with a 
low ionization parameter, the H II region has to be faint 
and old, and would therefore not have been observed 
in these surveys. For metallicities greater than about 
Z/Zq ~ 0.5, the models are highly degenerate in terms 
of both the metallicity and the cluster age. This accounts 



for the very tight distribution of observed points in the 
region of the plot where log([0 III]/H/3)< 0.0. 

Since the models of Dopita et al. (2000) were run, the 
EUV field predicted by the Starburst 99 v. 5 code has 
softened, thanks to improvements in the treatment of 
extended atmospheres. This leads both to the reduction 
in the limiting [N II] /Ha ratio seen in figure 3, which 
improves the fit with the observations, and the decrease 
in the limiting [O III]/H/3 ratio, which tends to make the 
fit somewhat worse at the low metallicity end. 

Note that the aging tracks make a sharp zig-zag caus- 
ing the 3 Myr isochrone to lie above the 2 Myr isochronc 
for metallicities above solar. This is particularly visible 
in panel (c) of figure 3, and is the result of the appear- 
ance of the Wolf-Rayet stars at ages ~ 3 Myr, which 
briefly harden the overall stellar radiation field before 
they themselves explode as Type II supernovae. 

The fit between theory and observation is much poorer 
when we consider the second of the Veilleux & Osterbock 
diagrams which plots the [S II] AA6717,3lA/H a ratio vs. 
the [O III] A5007A/H/3 ratio. This is shown in figure 
4. In general, the predicted [S II] lines are weaker than 
observed, by as much as a factor of two. A similar effect 
was seen in our earlier modeling, see figure 3 of Dopita 
et al. (2000), and is not fully understood. 

The [S II] AA6717,3lA/Ha ratio is more ambiguous 
than the [N II] A6584A/Ha ratio if we wanted to use it 
as an abundance indicator, particularly at high metal- 
licity. However, the 0.2 - 3 Myr isochrones still define a 
narrow range in line ratio space, albeit not exactly co- 
incident with the observations. The reasons for such a 
discrepancy between theory and observation in this, as 
well as in other line ratio plots will be discussed in section 
5. 

3.1.2. The BPT Diagram 

Baldwin, Phillips & Terlevich (1981) (here- 
inafter BPT) first proposed the use of the 
[O III] A5007A/[O II] AA3727,29A ratio vs the 
[O III] A5007A/H/3 line ratio as a diagnostic both 
to separate active nuclei from H II regions and also to 
provide an excitation sequence for H II regions. Our 
results for this line ratio are shown in figure 5. The 
[O HI]/[0 II] ratio is an excitation sensitive ratio, sensi- 
tive to the ionization parameter in the nebula. Because 
of their higher ionization parameters and their central 
stars with higher effective temperatures at a given age, 
the low-metallicity H II regions are characterized by 
systematically higher [O HI]/[0 II] ratios. Like all other 
ratios of optical forbidden lines to recombination lines, 
for a given stellar age, the [O III] A5007A/H/3 ratio 
starts off low at low metallicity, reaches a maximum 
at a certain metallicity, and then declines swiftly at 
higher metallicities. The reason for this was explained 
in the introduction in the context of the R23 ratio. As 
a consequence, the BPT diagram is highly degenerate 
in metallicity in the range 0.2 < Z/Zq < 2, explaining 
the tightness of the observed H II region sequence. This 
means, however, that the [O III]/H/3 ratio is of little use 
as an abundance diagnostic except when Z/Zq < 0.2. 
Even then it requires the observation of the [O III]/[0 II] 
ratio before an Oxygen abundance can be derived from 
it. 
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Fig. 3.— The [N II] A6584A/Ha cs. the [O III] A5007A/H/3 ratio. The abundances are coded by color; for Z/Zq0.05, 0.2, 0.4, 1.0 and 2.0 
the colors are brown, blue, green, black and red, respectively. In panel (a) we plot all models of a given metallicity, and in the remaining 
panels we plot ageing tracks at each metallicity and the corresponding isochrones for three likely values of the R parameter; log 1Z = 0, —2 
and -4, respectively. 



Although the theoretical sequence has much the same 
form as defined by the observations, it is notable that, 
for all values of the [O III]/[0 II] ratio, the maximum of 
the [O III]/H/3 lies below the maximum defined by the 
observations. Furthermore, even for the highest value of 
log 1Z, the zero-age models do not reach to the maximum 
[O III]/[0 II] ratio defined by the observations. As dis- 
cussed in section 5, below, both of these discrepancies 
can be ascribed to a single source, namely that the input 
spectrum we are using in the models is somewhat 'softer' 
than the EUV spectrum found in real H II regions. 



3.1.3. The R 23 Ratio 

As described in the introduction, the R23 ratio has had 
a long and checkered history as a metallicity indicator. 
However, the sensitivity of the [O III] A5007A/H/9 ratio 
to the cluster age shown in figure 4 suggests that we 
should plot the R23 ratio against the [O III]/[0 II] ratio 
in order to be able to investigate the age sensitivity of 
the R23 ratio. This diagnostic is shown in figure 5. 

For this diagnostic the fit between theory and observa- 
tion is quite poor, due mainly to the fact that the R23 ra- 
tio is systematically underestimated in the models. The 
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Fig. 4.— As figure 3 but for the second of the V&O diagnostics, the [S II] AA6717,3lA/Ha ratio vs. the [O III]/H/3 ratio. 



effect is larger than that seen in the [O III]/H/3 ratio (see 
figure 5) , showing that predicted strengths of the [O III] 
but especially the [O II] lines weaker than observations 
suggest. This will be discussed in the Discussion section, 
below. Points with low R23 would be explained as having 
metallicities greater than twice solar. 

This diagnostic plot shows immediately why the cali- 
bration of the R23 ratio in terms of metallicity has been 
so difficult. Not only is the range of R23 rather restricted, 
but is it degenerate in terms of metallicity over a rather 
generous range, 0.2 < Z/Z Q < 1.5, approximately. Not 
only that, but the ratio is fairly sensitive to excitation, 
as measured by the [O III]/[0 II] ratio. Because low 



metallicity H II regions are characterized by high ioniza- 
tion parameters, the low metallicity points tend to occur 
in the upper right of this plot, while, from the form of 
the isochrones at high metallicity, the spray of points to- 
wards the bottom left of this diagnostic can only be due 
to very high metallicity (Z/Zq > 2) H II regions with 
ages of 1-3 Myr. 

Based on these problems, we are constrained to con- 
clude that R23 is not a good abundance diagnostic. We 
will return to this point in Section 5, below. 

3.1.4. Dopita et al. (2000) Diagnostics 
Dopita et al. (2000) introduced two new abundance 
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Fig. 5.— As figure 3 but for the [O III] A5007A/[O II] AA3727.29A vs the [O III] A5007A/H/3 ratio. This is the diagnostic first used by 
Baldwin, Phillips & Terlevich (1981) (BPT) to separate AGN from H II regions. High mctallicity H II regions (Z/Zq > 2) are located in 
the lower left hand corner of this diagram. Note also how how the low metallicity models fall systematically below the main sequence of 
observed H II regions. 



diagnostics, based on the [N II] /[O II] line ratio. This 
ratio was favoured for two reasons which make it very 
sensitive to abundance. The first is that, over much of the 
range of interest, nitrogen is a secondary nucleosynthesis 
element, or at least has a large secondary component. 
Thus the ratio [N II] /[O II] will systematically increase 
with abundance. The second is connected with the fact 
that the [O II] AA3727,29A lines are found in the near UV 
part of the spectrum, whereas the [N II] A6584A is in the 
red part of the spectrum. This might be construed as a 
disadvantage in the use of the [N II] /[O II] ratio, since 
the reddening corrections are large, and correspondingly 



more uncertain. However, this problem is more than 
offset by the fact that the mean temperatures of H II 
regions are a strongly decreasing function of metallicity 
and at high abundances, the mean thermal energy of the 
electrons becomes too low to excite transitions with large 
electron excitation energy. Thus transitions such as the 
[O II] AA3727,29A lines are quenched, while the lower 
threshold energy [N II] A6584A continue to be excited. 
For this reason as well, the [N II] /[O II] line ratio falls 
with increasing abundance. 

The systematic decrease in the electron tempera- 
ture with increasing metallicity may also affect the 
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Fig. 6.— As figure 3 but for the R23— ratio vs. the [O III] A5007A/[O II] AA3727,29A ratio. In panels (b) to (d) the solar metallicity track 
has been omitted for clarity. This shows that the traditional R23 ratio is not a good abundance diagnostic in the range 0.2 < Z/Zq < 1.5, 
approximately, and that it is quite sensitive to the nebular excitation. The theoretical tracks are, however, a rather poor fit to the 
observations. This is because the models have [O III] AA4959,5007A and, even more so, [O II] AA3727,29A weaker than those observed. 



[O III] A5007A/[O II] AA3727,29A ratio. By itself, this 
would induce an increase in the [O III]/[0 II] ratio with 
metallicity. However, this is more than offset by two 
factors. Firstly, the mean effective temperature of the 
exciting stars falls with metallicity, and secondly, the 
mean ionization parameter is itself a decreasing function 
of metallicity. Both of these lead to a net decrease in the 
[O III]/[0 II] line ratio with metallicity. Nonetheless, the 
[O IIl]/[0 II] line ratio remains a sensitive indicator of the 
ionization parameter, which is strongly age dependent. 
Therefore, a plot of [N II] A6584A/[0 II] AA3727,29A 



against the [O III] A5007A/[O II] AA3727,29A ratio 
should separate stellar age from stellar abundance. This 
is illustrated in figure 7. 

It is clear that the [O III] A5007A/[O II] AA3727,29A 
ratio is an excellent diagnostic of both the stellar age 
and the nebular abundance. According to the theoretical 
tracks, all the observed H II regions have ages between 
and roughly 4 Myr, with most lying in the age range 
1-3 Myr. The division into abundance bins by use of the 
[N II] A6584A/[0 II] AA3727,29A is very clean. Provided 
that the calibration is correct, then the abundance can be 
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Fig. 7.— As figure 3 but for the [N II] A6584A/[0 II] AA3727.29A against the [O III] A5007A/[O II] AA3727,29A ratio used by Dopita et 
al. (2000). This is an excellent diagnostic of both metallicity and of stellar age. Note that the [O III] A5007A/[O II] AA3727,29A ratio is 
sensitive to both age and metallicity. The relationship between the observational points and the isochrones in this diagram suggests that 
most H II regions are observed in the age range 1-3 Myr. 



obtained using only the [N II] /[O II] ratio to a precision 
of roughly 0.15 dex. 

The [O III] A5007A/H/3 line ratio can be substituted 
for the [O III]/[0 II] ratio, since both are sensitive to ion- 
ization parameter and therefore to the age of the cluster. 
However, the separation into age classes is not so sharp 
using this ratio, as can be seen in figure 8. Furthermore, 
we have seen in figure 5 that the models tend to sys- 
tematically underestimate the [O III] A5007A/H/3 ratio, 
so therefore it follows that the age of the H II regions 
inferred using this diagnostic is also likely to be an un- 
derestimate of the true age. 



The third Dopita et al. (2000) diagnostic which utilizes 
the [N II] A6584A/[0 II] AA3727,29A ratio is that which 
plots this ratio against the [O III] A5007A/ [N II] A6584A 
ratio. The [O III] AA4959,5007A/ [N II] AA6548,84A ra- 
tio was first proposed by Alloin et al. (1979) as a more 
convenient abundance diagnostic, which is similar to the 
one we are using here. The observations and theoretical 
curves on this diagnostic are shown in figure 9. As can 
be seen in the diagram, both axes are quite sensitive to 
abundance, but the Alloin et al. (1979) ratio, by its use of 
the [O III] A5007A in the place of the [O II] AA3727,29A 
is also more sensitive to parameters other than the abun- 
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dance. 

Pettini & Pagcl (2004) have argued that the 
[O III] A5007A/ [N II] A6584A ratio is of itself a sen- 
sitive abundance indicator suitable to be used in the 
analysis of abundances in high-redshift galaxies. This 
is a little too simplistic. Because the [O III]/[N II] ratio 
is strongly dependent upon the excitation of the neb- 
ula, it is also sensitive to both the ionization parameter, 
as shown by Dopita et al. (2000) and to the age of the 
cluster of exciting stars. These sensitivities are as great 
as the sensitivity to abundance alone. Nonetheless, the 
[N II] A6584A/[0 II] AA3727,29A vs. [O III] A5007A/ 
[N II] A6584A diagnostic plot has great utility in sepa- 
rating abundance from cluster age, and the results that 
it gives are in broad agreement with figure 7. From these 
two plots we can conclude that most H II regions are ob- 
served in the age range 1-3 Myr with little sensitivity of 
this age range with abundance. 

This is in full accord with expectations. As shown in 
Paper II, after an age of 3 Myr, the absolute luminosity 
of H II regions in the recombination lines drops sharply, 
and due to the expansion of the H II region, the surface 
brightness drops even more rapidly. Both of these mili- 
tate strongly against the observation of old H II regions 
in extragalactic studies. Furthermore, very young H II 
regions tend to be dust obscured since they are still em- 
bedded in their placental molecular cloud complex. Such 
H II regions will not be bright at optical wavelengths. In 
Paper II we estimated the molecular cloud dissipation 
timescale for the molecular clouds to be about 1 Myr 
using the particular example of M17. The molecular 
cloud dissipation timescale was introduced in Paper I, 
and is the characteristic timescale r for the H II region 
to emerge from its compact or ultra-compact configura- 
tion, and is a measure of the solid angle (as measured 
at the exciting stars) which is intercepted by molecular 
clouds; f2(£)/47r = exp[— t/r). Clearly, the relative dearth 
of young H II regions beween 0.2 and 1.0 Myr, evident 
in figures 7 and 9, is consistent with r <~ 1 Myr. 

In conclusion, the [N II] A6584A/[0 II] AA3727,29A 
ratio vs. the [O III] A5007A/ [N II] A6584A provides a 
good diagnostic of both metallicity and age of the excit- 
ing cluster. The [O III] A5007A/ [N II] A6584A ratio can 
only be used by itself as an abundance indicator to the 
degree to which we can assume a (fairly restricted) char- 
acteristic age for the exciting clusters. Statistically, this 
may be valid, but the results of its use in single galaxies 
should be treated with caution. 

3.1.5. [N II] /Ha vs [S II] /Ha 

Because they lie very close to each other in wave- 
length, and because we have seen that [N II] is very 
sensitive to abundance, the [N II] A6584A/Ha vs. 
[S II] AA6717,3lA/Ha plot offers some hope as an abun- 
dance diagnostic. This has recently been investigated 
by Moustakas & Kennicutt (2005), who were interested 
to find out whether integrated spectra of galaxies were 
different from single H II regions when plotted in this 
way. We would expect to find such a difference in prin- 
ciple, because the integrated spectra do not discriminate 
against the emission from old H II regions in the man- 
ner described in the previous section. The results of our 
modelling is shown in figure 10. 



In Paper II, we found that the diffuse component due 
to old, evolved, large and faint H II regions could amount 
to 20-30% of the total Ha emission, and suggested that 
this could be identified with the Diffuse Ionized Medium 
(DIM) which is almost ubiquitously seen in disk galax- 
ies (Hoopes, Waterbos & Greenwalt 1996; Martin 1997; 
Wang, Heckman & Lehnert 1998, 1999). This is usually 
refereed to as the "Reynolds Layer" in the context of our 
own Galaxy; see a recent review by Reynolds (2004). In 
the past, the excitation of the DIM has been usually as- 
cribed to leakage of ionizing photons from young H II 
regions. However, the emission from large, faint, low 
pressure and evolved H II regions might provide a better 
explanation. 

In figure 10, the models of old H II regions with 
mctallicities Z/Zq > 1.0 reach up to line ratios of 
[S II] AA6717,3lA/Ha - 1 and [N II] AA6548,84A/Ha 
~ 2/3, which is certainly consistent with observed line 
ratios encountered in the DIM. However, in our Galaxy, 
line ratios even greater than these are seen. To explain 
these, we would have to invoke an additional heating 
mechanism, such as shocks, or a harder EUV field pos- 
sibly produced by old ex-nuclei of Planetary Nebulae, 
which are sufficiently numerous to produce an apprecia- 
ble EUV field out of the plane of the Galaxy. 

3.2. Infrared Line Diagnostics 

The ISO satellite first demonstrated the utility of mid- 
IR line diagnostics in the analysis of H II region spectra. 
For normal galactic H II regions, Giveon et al. (2002) 
obtained extensive ISO spectroscopy, as did Peeters et 
al. (2002) in the context of Ultra-Compact (UC) H II 
regions. Later, Verma et al. (2003) measured integrated 
spectra for a number of famous Starburst galaxies. We 
will use these data here, although we note that both the 
high-resolution and low-resolution modes of the Infrared 
Spectrometer (IRS) (Houck et al. 2004) on the Spitzer 
Space Observatory are producing much new and excel- 
lent data on both individual H II regions as well as on 
spectroscopy of whole galaxies. 

Within the spectral range of the IRS, such line ra- 
tios that can be constructed are primarily sensitive to 
the nebular excitation. Examples include the [Nc III] 
15.5yum / [Ne II] 12.8/im ratio, the [S IV] 10.5^m / [S III] 
18.7/im ratio, the [S IV] 10.5/xm / [Ar III] 9.0^m ratio, 
and the [Nc III] 15.5^m / [S III] 18.7^m ratio. Unfor- 
tunately, these arc of little use in measuring anything 
other than the hardness of the EUV radiation field. All 
of the strong lines accessible in this region of the spec- 
trum are produced by a— process elements, and there- 
fore the abundance ratios stay approximately constant 
with abundance. All the lines have low energy thresholds 
for excitation, and therefore none of them are particu- 
larly sensitive to the electron temperature, which varies 
strongly with abundance. In figure 11 we present only 
one example to demonstrate how poor these line ratios 
are as abundance diagnostics when used alone. 

Much of the data used in this plot refers to observations 
of ultra-compact H II regions, which are not exactly com- 
parable to "normal" H II regions. In particular, they are 
extremely dense, most are excited by only a single OB 
star, and many of the ionizing photons from the exit- 
ing star are absorbed by dust rather than by the ionized 
plasma. This accounts for much of the offset between 
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Fig. 8.— As figure 3 but for the [N II] A6584A/[0 II] AA3727,29A against the [O III] A5007A/H/3 ratio, as used by Dopita et al. (2000). 
This is also a good diagnostic of both metallicity and of stellar age. Given that the models tend to underestimate the [O III] A5007A/H/3 
ratio, the age of the H II regions using this diagnostic is also likely to be an underestimate. 



the (low density) models and the data in figure 11. This 
effect is particularly important for the [S IV] 10.5/im / 
[S III] 18.7/im ratio. These points, and this particular 
diagnostic plot was thoughly discussed in our paper on 
compact H II regions (Dopita et al. 2006b). 

In order to form a line ratio that is sensitive to abun- 
dance, we need to use a ratio with a Hydrogen or He- 
lium recombination line. Unfortunately, these are rather 
faint in the mid-IR. About the best line we have is 
Bra A4.051/im. Alternatively, Pfund— a A7.458/im might 
be useable, although it is often difficult to observe against 
the PAH features seen at this wavelength. In figure 12 we 



show one such diagnostic, the [Ne III] 15.5/im / [Ne II] 
12.8/im ratio versus the [Ne III] 15.5/im / Bra A4.051/im 
line ratio. This has the advantage that, at a given value 
of the excitation measured by the [Ne III] 15.5/im / 
[Ne II] 12.8/xm ratio, the [Ne III] 15.5/im / Bra A4.051/mi 
line ratio is a monatonic function of abundance up to at 
least Z ~ 2.0 x Z d t, and so it provides a clean abun- 
dance diagnostic, unlike the comparable optical line ratio 
shown in figure 5. 

Currently, the uncertainties in the observational data 
are to large for this to provide a very accurate abundance 
diagnostic. However, we may note that the observations 
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Fig. 9.— As figure 3 but for the [N II] A6584A/[0 II] AA3727,29A against the abundance-sensitive [O III] A5007A/ [N II] A6584A ratio, 
as proposed by Alloin et al. (1979). Although the [O III]/[N II] ratio is indeed sensitive to abundance, it is also very sensitive to ionization 
parameter as shown by Dopita et al. (2000) and therefore to stellar age. This diagram suggests that most H II regions are observed in the 
age range 1—3 Myr, in good agreement with age estimate obtained by figure 7. 



from Giveon et al. (2002) are consistent with abundances 
between 2.0 and 0.4 Zq, and with ages between and 
2 Myr. This abundance range is exactly what we would 
expect to find for Galactic H II regions over the range 
of Galactocentric radii observed by him. It is interest- 
ing that the ages are lower than those derived for H II 
regions at optical wavelengths. This is consistent with 
the selection criteria for IR-bright, compact H II regions. 
This selects for young H II regions. Additionally, the use 
of mid-IR wavelengths allow the H II region emission to 
be seen almost unimpeded by the surrounding molecular 



material, so the molecular cloud dissipation timescale is 
not a relevant parameter here. 

4. INTEGRATED SPECTRA 

The diagnostics presented in the previous sections were 
designed to derive parameters for single H II regions. 
However, what we seek to do now is to develop a spec- 
tral synthesis by taking recombination line flux averaged 
spectra along an evolutionary track of an H II region with 
given parameters. This reduces the number of free pa- 
rameters to two; Z/Zq and log 1Z. The observations can 
then be used to constrain the values log??., allowing the 
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Fig. 10.— As figure 3 but for the [N II] A6584A/Ha vs. [S II] AA6717,3lA/Ha diagnostic line ratios. As in Figure 3 the data for 
individual H II regions are shown as crosses, and the integrated spectra from Moustakas & Kennicutt (2005) are shown as points. Because 
these are both ratios of lines in the red portion of the spectrum, and because N is partially a secondary element, this diagnostic turns out 
to be a good abundance diagnostic for Z/Zq < 1.0. The models of the older and more diffuse H II regions are located close to the right 
hand boundary of the observational points. This gas can be identified with the warm ionized medium (WIM) in galaxies. 



extragalactic abundance sequence for ensemble averages 
of evolving H II regions to be derived. This technique 
opens the way to the accurate study of the metallicity 
evolution of both individual galaxies and the Universe. 

The model we adopt assumes that clusters are formed 
continuously, and that any single cluster has all of its 
stars born coevally. As described above, the cluster mass 
function does not enter explicitly into the models, which 
are characterized by their R— values. The process of pro- 
ducing an integrated spectrum for any Z is therefore 



simply to integrate the fluxes from a single H II model 
through time, until it has faded such that > 96% of all 
ionizing photons have been emitted. This limits the inte- 
gration time to lie within the first 5 — 6.5 Myr. Stochastic 
effects due finite sampling of the stellar IMF - important 
at low cluster masses - are not included. However, the 
form of the cluster initial mass function either as inferred 
in Paper II, or derived from stellar number counts is such 
that each logarithmic bin of cluster mass contributes a 
similar amount to the total flux. Massive clusters will 
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Fig. 11.— The [Ne III] 15.5/mi / [Ne II] 12.8/xm ratio versus the 
[S IV] 10.5/mi / [S III] 18.7/mi ratio. H II regions from Giveon 
et al. (2002) are plotted as filled circles, UC H II from Peeters et 
al. (2002) as open circles, and Starburst galaxies from Verma et 
al. (2003) as filled squares with error bars. The theoretical curves 
create as superb an example of theoretical spaghetti as anything 
these authors have seen. 
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Fig. 12.— The [Ne III] 15.5/xm / [Ne II] 12.8/im ratio versus 
the [Ne III] 15.5/xm / Bra A4.051/xm line ratio. H II regions from 
Giveon et al. (2002) are plotted as filled circles. This diagnostic 
provides a clean abundance determination which is monotonic over 
the full abundance range. The isochrones are given for the case 
log Ti = —2. This diagram is an infrared analog of fig 5, but is 
a much cleaner diagnostic because these IR line ratios are very 
insensitive to the nebular temperature. 



therefore tend to wash out the stochastic effects of low- 
mass clusters. 

On line ratio diagnostics, there is no simple relation- 
ship between the ensemble averaged line ratio and the 



line ratios in the individual H II regions which make up 
this average. For example, a high ionization line such 
as [O III] A5007A will be mostly produced in the very 
youngest H II regions with the highest R— values. Low 
ionization lines such as the [S II] AA6717,3lA lines will 
have a much greater weighting from older evolved H II 
regions with low R— values. Indeed, for this line the ma- 
jority of the line flux may come from regions lying outside 
the bright H II regions that would be normally selected 
for observation in studies of individual H II regions - in 
other words, from regions that form part of the 'warm 
ionized medium' (WIM). 

Given these quite different selection biases, it is in- 
teresting to compare line ratio observations of the in- 
tegrated line fluxes of galaxies with those of individual 
H II regions. Such a comparison helps us to quantify the 
importance of the WIM in determining the global line 
ratio. For this purpose, we have used the observations 
of individual H II regions listed and used above ( 3 and 
references in Section 3.1), and we have drawn our sam- 
ple of integrated line fluxes of galaxies from the Nearby 
Field Galaxy Survey (NFGS) (Jansen et al. 2000), and 
the recent survey by Moustakas & Kennicutt (2005). 

In figure 13 we plot the main V&O diagnostic 
[N II] A6584A/Ha vs. [O III] A5007A/H/3. This should 
be compared with figure 3. The models provide a good 
fit to the integrated spectra of galaxies. Indeeed, the fit is 
better than that of the models of individual H II regions 
compared to data in figure 3. This is especially appar- 
ent at the low abundance end (where most observations 
are of relatively nearby dwarf irregular galaxies such as 
NGC6822, IC1613, IIZw40 and the like). This suggests 
that the observations of these low-abundance H II re- 
gions may be affected by aperture effects. This would be 
caused by observers picking out bright high-excitation 
sub-regions of the H II region from their long-slit spectra 
for analysis, rather than observing across the whole H II 
region to obtain an integrated spectrum of it. This would 
preferentially enhance the [O III] A5007A/H/3 ratio while 
at the same time depressing the [N II] A6584A/Ha ratio, 
moving the observed point in line ratio space towards the 
top and the left. 

The fit with the second V&O diagnostic, the 
[S II] AA6717,3lA/Ha vs. [O III] A5007A/H/3 (figure 14) 
is less convincing. The integrated spectra are displaced 
to the right on the [S II] AA6717,3lA/Ha ratio compared 
with the spectra from individual H II regions, showing 
the importance of the contribution of the more diffuse 
ionized component. It is interesting to note that most of 
the highest points in terms of the [O III] A5007A/H/3 ra- 
tio are contributed by the Moustakas & Kennicutt (2005) 
data set. In this respect there are systematic differences 
between this data set and the NFGS survey. The reason 
for this is not apparent to us. 

In the models, the [N II] AA6548,84A are the ma- 
jor coolants in the outer parts of the H II regions. It 
is possible that the abundance of Nitrogen in the ion- 
ized gas is somewhat overestimated (by an underesti- 
mate of the depletion, for example), leading to saturation 
in the [N II] A6584A/Ha ratio, and a reduction in the 
[S II] AA6717,3lA/Ha ratio. The problem with the mod- 
els in reproducing the observed [S II] AA6717,3lA/Ha 
ratios is systematic. Even with the old Starburst 99 v2 
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models, a similar problem was encountered by Dopita 
et al. (2000). This will be discussed further in the next 
section. 

Let us now turn our attention to the BPT diagram, 
the [O III] A5007A/[O II] AA3727,29A ratio vs. the 
[O III] A5007A/H/3 ratio (Baldwin, Phillips & Terlevich 
1981). This is plotted in figure 15. Once again, the fit 
between the theory and observation is good, although on 
this diagnostic the data once again fall somewhat high 
relative to the theory in the [O III] A5007A/H/3 ratio at 
the low abundance end Z/Z & < 0.4. 

We note that the [N II] A6584A/Ha ratio covers 
quite a large range and displays considerable sensitiv- 
ity to abundance. However, the same is not true of the 
[O III] A5007A/H/3 ratio. Only two points in the NFGS 
survey lie below log[N II] A6584A/Ha = -1.5, and since 
the Z/Z Q curve lies to the left of log[N II] A6584A/Ha 
= —1.5, figure 13 suggests that the lower abundance 
bound to the NFGS survey is about Z/Z Q < 0.2. This is 
also consistent with what would be implied from figure 
15. Both figure 14 and figure 15 suggest that the upper 
bound is about Z/Z & ~ 2. 

The first of the Dopita et al. (2000) diagnos- 
tics is shown in figure 16. This plots the 
[N II] A6584A/[0 II] AA3727,29A ratio against the 
[O III] A5007A/[O II] AA3727,29A ratio, and should be 
compared with figure 7. There is very little distinction 
between the distribution of the integrated spectra from 
the NFGS survey and the observations of individual H II 
regions on this diagram. Essentially any value of \ogTZ 
seems to be permitted, although the majority of the 
points are consistent with a range to -4 in this parame- 
ter. This is consistent with what we might have expected, 
given that the parameters which enter into R , M c \ and 
P/k, cover at least the range 100 < M ci /M Q < 10 5 and 
10 4 < P/k < 10 7 cur 3 K. 

The abundance range of the NFGS survey suggested 
by this plot is broader than that given by the V&O or 
BPT diagnostics, since the theory implies the existence 
of some galaxies with very low abundances. However, 
we must beware once again of the observational prob- 
lems which enter into plots like this. The [N II] A6584A 
line and the [O II] AA3727,29A lines lie at opposite ends 
of the spectrum, and so there is a very large redden- 
ing correction applied to the data. Generally speaking a 
standard reddening law is used to do this. In the case 
of the NFGS survey, for example the Clayton, Cardelli 
& Mathis (1989) extinction curve was used. However, 
the foreground dust screen in front of an individual H II 
region (let alone a whole galaxy) has a complex fractal 
structure, and rather than using an extinction curve, an 
attenuation curve should be used instead, see Fischera, 
Dopita & Sutherland (2003); Fischera & Dopita (2005) 
and Paper I. 

The effect of an attenuation law will be to cause the 
strengths of lines in the red to be underestimated, while 
UV lines such as the [O II] AA3727,29A lines will be sys- 
tematically overestimated. These effects will mean that 
the observations corrected using an extinction law will 
land systematically too far to the left and too low on this 
diagram, while on (almost) extinction-free diagnostics 
such as the V&O plots, the data will plot onto the correct 
position. If this effect is the cause of the displacement 
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Fig. 13. — The synthesised H II region spectra on the V&O di- 
agnostic [N II] A6584A/Ha vs. [O III] A5007A/H/3 (c.f. figure 3). 
The integrated galaxy spectra from the Nearby Field Galaxy Sur- 
vey (NFGS) (Jansen et al. 2000) arc shown as filled circles, while 
the crosses are for individual H II regions drawn from the sources 
referenced above, to show the degree to which the integrated spec- 
tra differ from the single source spectra. Note the large disparity 
between the integrated spectra and the observations of individual 
H II regions at the low abundance end (Z/Zq < 0.4). Compared 
with earlier models, the theory provides a good fit to the inte- 
grated galaxy data providing that log 1Z > —4, and emphasizes the 
degeneracy of this V&O diagnostic. 
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Fig. 14. — The synthesised H II region spectra on the V&O di- 
agnostic [S II] AA6717,3lA/Ha vs. [O III] A5007A/H/3 (c.f. figure 
5). The integrated galaxy spectra are shown as filled circles, and 
include both the NFGS data and the data from Moustakas & Kcn- 
nicutt (2005). Both the [S II] AA6717,3lA and the [O III] A5007A 
lines seem to be a little weak in the models compared with the 
observations. 
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Fig. 15. — As figure 13, but for the BPT diagram, the 
[O III] A5007A/[O II] AA3727,29A ratio vs [O III] A5007A/H/3 
. This should be compared with figure 5. 



between theory and observation seen in figure 15, this 
suggests that the [O III] A5007A/[O II] AA3727,29A ratio 
may be systematically underestimated by ~ 0.3dex. This 
would imply that the [N II] A6584A/[0 II] AA3727,29A 
ratio may be underestimated by a similar amount, since 
the extinction is determined from the Ha/H/3 ratio which 
means that the observed [N II] A6584A/[0 III] A5007A 
ratio must necessarily be correct. The effect of a correc- 
tion of this magnitude is shown on figure 16 as an arrow 
marked 'extinction'. With a correction of this magni- 
tude, not only can the lower bound of the abundance in 
the NFGS survey be reconciled with that given by figure 
13 (Z/Zq ~ 0.2), but also with the upper bound indi- 
cated by figure 15 (Z/Z@ ~ 2). The inferred range on 
log 1Z would then be 2 to -4. 

A systematic overestimate of the [O II] AA3727,29A 
lines would also lead to a systematic overestimate of R23. 
This may well account for much of the discrepancy be- 
tween theory and observation in figure 6. 

A similar effect is seen in the second of the Dopita et 
al. (2000) diagnostics is shown in figure 17. This plots 
the [N II] A6584A/[0 II] AA3727,29A ratio against the 
[O III] A5007A/H/3 ratio. Once again, a systematic ex- 
tinction correction error of 0.3 dex would shift the obser- 
vational points into a range which is consistent with the 
both the abundance and R bounds inferred above. 

Finally, we plot the [S II] AA6717,3lA/Ha ratio versus 
the [N II] A6584A/Ha line ratio in figure 18 (compare 
with figure 10). This shows, once again, the value of the 
[N II] A6584A/Ha line ratio as an abundance diagnos- 
tic in its own right. It also illustrates, once again, our 
difficulties with the [S II] AA6717,3lA/Ha ratio, with 
the observations, particularly the Moustakas & Kenni- 
cutt (2005) data set falling to the right of the theoretical 
curves in this ratio. 

Note how the integrated spectra on figure 18 are con- 
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Fig. 16. — As figure 13 but for the 

[N II] A6584A/[0 II] AA3727.29A ratio against the 
[O III] A5007A/[O II] AA3727,29A ratio used by Dopita ct 
al. (2000). This should be compared with figure 7. This is perhaps 
the most unambiguous diagnostic for the determination of both 
mctallicity and log7£. For the meaning of the arrow marked 
'extinction' see the text. 
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Fig. 17.— As figure 13, but for the 

[N II] A6584A/[0 II] AA3727.29A ratio against the 
[O III] A5007A/H/3 ratio used by Dopita ct al. (2000). This should 
be compared with figure 8. This is a good diagnostic for the 
determination of both mctallicity and log TZ. The observations and 
the theory would give consistent answers in all integral spectrum 
diagnostics if a systematic correction is made to the observations 
of the magnitude illustrated by the arrow marked 'extinction' in 
this and the previous figure (see text). 
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Fig. 18.— As figure 13, but for the [N II] A6584A/Ha ratio 
against the [S II] AA6717,3lA/Ho ratio. This should be compared 
with figure 10. This once again illustrates our problems in mod- 
elling the [S II] line strengths. 



centrated along the right-hand side of the region defined 
by the observations of single H II regions. This reflects 
the fact that integrated galaxy spectra include the dif- 
fuse component of the ISM, even when this has a low 
surface brightness, because the total flux in the diffuse 
component may be an appreciable portion of the total 
emission. Indeed, a recent study by Oey et al. (2006) 
suggests that the "diffuse" fraction may exceed 50% in 
nearby starburst galaxies. 

5. DISCUSSION & CONCLUSIONS 

It should now be clear that these models represent a 
significant advance in our understanding of the strong 
line spectra of H II regions. In particular, the models ex- 
plicitly account for the dynamical evolution of the H II 
region controlled through the mechanical energy input 
by both the stellar winds and supernova explosions. Two 
factors largely control the excitation of the young H II 
regions, both of these dependent on the chemical abun- 
dance. These are: 

• The effective temperature of the exciting stars, 
which increases with decreasing metallicity, and 

• The control of the ionization parameter by the stel- 
lar winds in the expanding H II regions. This 
was shown in Paper II to increase with decreasing 
metallicity 

Since both of these effects operate to increase the excita- 
tion of low abundance H II regions, we therefore find an 
abundance dependent excitation, which much improves 
the fit of the theory compared to models in which the 
ionization parameter is treated as a free variable. In our 
models the ionization parameter is replaced by the R 
parameter, the ratio of the cluster mass in solar masses 



to the pressure in the ISM; P/k (cm -3 K). This param- 
eter is used because all models with a given R have a 
unique relationship between their instantaneous ioniza- 
tion parameter and time. This allows us to construct 
evolutionary tracks for ensembles of H II regions. 

With increasing age, the ionization parameter falls al- 
most monotonically, and the H II region fades after reach- 
ing its maximum brightness in Ha after 1-2 Myr. The 
presence of Wolf-Rayet stars after about 3.5 Myr causes 
the effective temperature of the cluster to briefly increase. 
This increases the excitation for a short time, causing 
the 3 Myr and the 4 Myr isochrones to cross over each 
other on diagnostic plots for the cases of high metallicity 
{Z/Zq > 0.5, approximately). 

In this paper, we have identified a number of useful di- 
agnostics for deriving the abundance and either the age, 
in the case of isolated H II regions, or the R parame- 
ter, in the case of integrated spectra of galaxies. Our 
diagnostics for the integrated spectra are the first which 
explicitly take into account the fact that when we observe 
a whole galaxy, we are observing an ensemble average of 
H II regions of all ages, sizes and central cluster masses. 
The simplification we have been able to make by using 
the R parameter has rendered this problem tractable to 
a more simple analysis other than treating each H II re- 
gion separately, and adding the whole lot together at the 
end. 

In comparing the models to the observational data we 
have identified some potential problems with the way 
that observations of H II regions have been made and 
analyzed. In particular: 

• Long-slit observations do not integrate over the 
whole of an H II region, particularly in the case 
of observations of nearby, low metallicity dwarf ir- 
regular galaxies. By picking out the brightest re- 
gions of such H II regions, we bias the data towards 
the high-excitation regions and the spectrum is not 
representative of the whole H II region (Luridiana 
& Peimbert 2001; Stasinska & Izotov 2003) and 

• The use of a reddening correction using an extinc- 
tion curve, rather than an attenuation curve (which 
is appropriate for spatially extended objects) will 
cause the line intensities of UV lines to be system- 
atically overestimated, biasing the strong line di- 
agnostics which use the [O II] AA3727,29A A lines. 

This second effect would reveal itself in the spectro- 
photometry as a failure to reproduce the correct intensi- 
ties of the higher members of the Balmer Series, such 
as or He, and so should be quite easy to check 
for in observational data sets. Such an effect would 
also influence the measurement of the strength of the 
[O III] A4363A line, again causing it to be overestimated. 
This could be part of the reason for the long-standing 
discrepancy between the abundances derived from the 
strong lines and those obtained from measurements of 
the electron temperature using ratios such as [O III] 
A4363A/[0 III] A5007A , see Bresolin, Garnett & Kcnni- 
cutt (2004) for a recent discussion. We expect that the 
major part of the discrepancy will turn out to be caused 
by temperature fluctuations. However, this could be eas- 
ily corrected for by assuming the Case B recombination 
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line ratio between H/3 and H7, and measuring the forbid- 
den line strengths relative to the nearest Balmcr lines in 
the spectrum. 

There is clearly a need for integral field spectropho- 
tometry on extragalactic H II regions. This would allow 
efficient construction of integrated spectra and assist in 
the removal of stellar flux from the spectrum. Such ob- 
servations would not only help to remove the remaining 
discrepancies between theory and observation but also 
give insight into the use of theoretical integral spectra in 
the analysis of the chemical evolution history of the early 
Universe. 

Our models have their deficiencies. In general, 
all the low ionization species are predicted weaker 
than observed. These include the [O I] A6300A, the 
[S II] AA6717,3lA, and the [O II] AA3727,29A, although 
these lines are probably affected by the reddening cor- 
rections made to the observations. 

What could enhance the [S II] AA6717,3lA/Ha ratio 
or the [O I] A6300A/Ha ratio? First, both the the O I 
and the S II ion are very much more sensitive than the 
N II ion to the diffuse radiation field in H II regions. 
They are therefore much enhanced in the vicinity of ion- 
ization fronts, and in shadow regions behind elephant 
trunk features, which are obscured from direct stellar 
radiation and illuminated only by the diffuse nebular ra- 
diation field. Second, both are much enhanced in shock 
regions. Indeed high [S II] /Ha ratios are used as a means 
of identifying supernova remnants cither embedded in 
H II regions or in the ISM at large, and [O I] /Ha ratios 
reach 0.5 or higher in shock-excited regions. Neither of 
these effects are taken into account in the models. 

In addition, the O I ion is locked by charge-exchange 
reactions to the concentration of H I. Furthermore, the 
collision strength of [O I] A6300A increases with electron 
temperature. To produce a high [O I] A6300A/Ha ratio 
in a purely photoionized region requires the presence of 
a hard radiation field. This could originate either as a 
thermal soft X-ray continuum from the shocked bubble of 
stellar wind material, or it could simply be an indication 
of a hotter sub-component in the global stellar radiation 
field. 

If the radiation field produced by Starburst 99 v5. is 
too soft, this could be rectified in at least three ways, 
First, by the use of a flatter IMF in the upper mass ranges 
(M/Modot > 10). This could be possible because obser- 
vationally we have very few constraints other than by use 
of the excitation of H II regions. Second, by adding an- 
other stellar component such as mass-exchange massive 
binaries to the stellar synthesis. 

A third - and quite likely possibility is that the EUV 
spectrum of the massive stars used in Starburst 99 v5. 
may be in error. With previous versions of Starburst 
99 based on the CoStar models (Schaerer & de Kotcr 
1997), the EUV spectrum was, rather, too hard. How- 
ever, Morisset et al. (2004) has demonstrated that the 
theoretical number of ionizing photons emitted by the 
central star and the shape of the ionizing spectrum is 
highly dependent upon the atmospheric model used. At 
one extreme are the CoStar models, at the the other 
extreme are found the 'classical' Kurucz (Kurucz 1991, 
1994) plane-parallel LTE line blanketed models used in 
Starburts99 v5.. The Kurucz models display the softest 



EUV spectra and the lowest photoionising flux, and are 
restricted to log g = 3.0 for the models with higher ef- 
fective temperature, T c g. In between these extremes are 
found three sets of models with fairly similar photoinis- 
ing flux predictions (Morisset et al. 2004), the TLUSTY 
models of Lanz & Hubeny (2003a, b), the WM-Basic 
models of Pauldrach, Hoffman, & Lennon (2001) and the 
CMFGEN models Hillier & Miller (1998). Of these, the 
TLUSTY models are plane-parallel hydrostatic, while 
the other two arc spherical dynamic atmospheres. All 
of this offers plenty of potential for error in the estima- 
tion of the the hardness of the EUV spectrum for a given 
cluster mass and metallicity. 

Given that the [O III] A5007A/H/3 ratio of the mod- 
els also appears to be on the low side, we believe that 
an increase in the stellar effective temperature, either by 
a flattening of the IMF or by some other means offers 
the best chance to more properly match the theory with 
the observations. The models are too weak in their for- 
bidden line strengths overall. The Stoy (1933) method 
(Pottasch & Preite-Martinez 1983) of deriving the effec- 
tive temperature of the exciting stars relies upon mea- 
suring the total flux cooling lines (which is a measure of 
the mean energy per photoionization) to a recombina- 
tion line (which counts the number of photoionizations). 
Thus, if we want to increase the absolute strengths of the 
forbidden lines, we have to increase the stellar effective 
temperature. This increase would also strongly increase 
the [O I] A6300A/Ha ratio by both increasing the ex- 
tent of the partially-ionized zonr and by increasing the 
electron temperature within it. 

A niggling uncertainty with the modeling is the mis- 
match between the Starburst 99 v5 abundance set and 
the abundance set given in table 1. The Starburst 99 v5 
has to rely on the older stellar evolution models which 
use the old value of the solar abundances. In an at- 
tempt to quantify this effect, we ran a test model with a 
OAZq spectral synthesis cluster model from Starburst 99 
v5 embedded in a nebula with I.OZq (using the abun- 
dances given in table 1). This made surprisingly little 
difference, changing critical line ratios by 0.1 dex or less, 
except for the [O I] A6300A/Ha ratio, which increased by 
34%. We therefore conclude that the results presented 
in this paper are reasonably secure against any future 
adjustment of stellar atmospheric abundances. 

In conclusion, the major advances in modelling made 
in this paper are: 

• For individual HII regions, it becomes possible to 
estimate the ages of the exciting stars from the po- 
sitions of the observations on HII region isochrones. 
Previously, this was not possible, since the ioniza- 
tion parameter and the effective temperature of the 
cluster (determined by the ageing of its stars) pro- 
duce similar effects on the emitted spectrum. The 
self-consistent treatment we have developed which 
includes the dynamical evolution of the HII region 
consistent with the properties of the stellar wind 
generated by the central cluster, and which recog- 
nises that the stellar wind determines the pressure 
in the HII region allows us to eliminate the use of 
arbitrary geometries. 

• In addition, this work gives a methodology to con- 
struct much more "realistic" integrated spectra of 
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Table 2. Line ID and Wavelength list 



Line ID # Ion A (A) Line ID # Ion A (A) 



01 


HI 


1215.7 


02 


CIII] 


1909.1 


03 


CIII] 


1911.2 


04 


[C II] 


2325.2 


05 


Mgll 


2797.9 


06 


[O II] 


3726.0 


07 


[O II] 


3728.7 


08 


HI 


3797.9 


09 


HI 


3835.4 


10 


[Ne III] 


3868.7 


11 


He I 


3888.6 


12 


HI 


3889.1 


13 


[Ne III] 


3967.4 


14 


HI 


3970.1 


15 


H I 


4104.7 


16 


H I 


4340.5 


17 


He I 


4471.5 


18 


HI 


4861.3 


19 


[O III] 


4958.8 


20 


[O III] 


5006.8 


21 


He I 


5875.6 


22 


[O I] 


6300.2 


23 


[N II] 


6548.0 


24 


H I 


6562.8 


25 


[N II] 


6583.3 


26 


He I 


6678.2 



27 


[S II] 


6716.3 


28 


[S II] 


6730.7 


29 


[Ar III] 


7135.7 


30 


[Ar III] 


7751.0 


31 


[S III] 


9069.3 


32 


H I 


10049.4 


33 


HI 


10049.4 


34 


He I 


10830.0 


35 


He I 


10833.0 


36 


HI 


10938.1 


37 


HI 


12818.1 


38 


HI 


18751.0 


39 


HI 


26252.0 


40 


HI 


40512.0 


41 


[Ar II] 


69832.8 


42 


[Ar III] 


89892.5 


43 


[S IV] 


105221 


41 


[Ne II] 


128115 


45 


[Ne III] 


155513 


46 


[S III] 


186821 


47 


[S III] 


336366 


48 


[Si II] 


347941 


49 


[O III] 


517972 


50 


[N III] 


573845 


51 


[O III] 


883017 



HII galaxies, something that has not been possible 
before in the approximations which assume that a 
whole galaxy can be represented by a single spher- 
ical H II region with a chosen abundance and ion- 
ization parameter. 
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APPENDIX 

Here we give a summary of the strong line flux ratios with respect to H/3 as computed in the photoionization models. 
The line identifications, the ionic species and the wavelengths of each line listed in the electronic version of this paper 
are given in Table 2. 

In the electronic version of the paper, separate tables are given for each of the abundances computed; 0.05, 0.2, 0.4, 
1.0 and 2.0 Zq, and in each of these tables, line fluxes with respect to H/3 = 1.0 are given for each of the 51 lines 
listed in Table 2. These spectral line intensities are given for each value of log TZ computed; +2, 0, -2, -4 and -6. The 
quantity R is that defined in the text of the paper: 



K = 



M G 



P/k 



10 4 cm" 3 K 



(1) 



The electronic tables (Tables 3 through 7) will be available when the paper is published. 
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